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Abstract. We start from an interpretation of the _BC2 -symmetric "Type I" 
(elliptic Dixon) elliptic liypergeometric integral evaluation as a formula for a 
Casoratian of the elliptic hypergeometric equation, and give an extension to 
higher-dimensional integrals and higher-order hypergeometric functions. This 
allows us to prove the corresponding elliptic beta integral and transformation 
formula in a new way, by proving both sides satisfy the same difference equa- 
tions, and that the difference equations satisfy a Galois-theoretical condition 
ensuring uniqueness of simultaneous solution. 



1. Introduction 

Plain hypergeometric functions and their g-analogues are widely used in math- 
ematics and mathematical physics. They can be defined either as infinite series or 
contour integrals, more or less on an equal footing As to the recently discov- 
ered elliptic hypergeometric functions, the situation with them is different — their 
general instances are defined only via integral representations. The general concept 
of elliptic hypergeometric integrals was introduced in [T4j. Elliptic beta integrals 
[H [ini Uni [m are the simplest representatives of integrals of such type. In the 
univariate setting there is only one elliptic beta integral JT3] , presently the top level 
known generalization of the Euler beta integral. In the multivariable case such 
integrals are grouped in three classes. 

The n-dimensional type I elliptic beta integrals contain 2n + 3 free parameters, 
and there are two known general methods of proving them [10|,ll6j. The type II inte- 
grals have a smaller number of parameters, and admit a straightforward derivation 
from the type I integrals f?, 13^ . Both types of these exactly computable integrals 
admit higher-order extensions with more parameters, with associated transforma- 
tion laws [To]. For the BCn root system, these are elliptic analogues of an integral 
due to Dixon [5] (type I) and of Selberg's famous integral [2] (type II, with 5-1-1 
parameters), respectively. Multiple elliptic beta integrals of the third class [Hj can 
be represented as determinants of univariate integrals, which, in turn, reduce to 
computable theta function determinants. 

In the present paper, we follow up on the observation implicit in [lOj that the 
elliptic Dixon integrals can be expressed as determinants of univariate elliptic hy- 
pergeometric integrals, higher-order analogues of the elliptic beta integral. This 
allows us to give a new proof of the corresponding evaluation formula and of the 
related transformation formula established first in [10] . The evaluation result lifts 
Varchenko's determinant of univariate plain hypergeometric integrals [19] and the 
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Aomoto-Ito determinant [3] to the elhptic level and enriches the list of known com- 
putable determinants compiled in [7]. 

We use the following notation. The key infinite product is 

oo 
k=0 

where \p\ < 1 and z G C. The elliptic theta function has the form: 

ep{z) := {z-p)oo(z^^;p)co, 
where z G C*. It obeys the properties 

Opipz) = ^p(z-i) = -z-'ep{z) 
and 9p{z) — for z ~ p^. We follow the standard useful convention that 

m 

9p{ai, ...,a^):=Y[ ^p(«fc)' - Op{tz)ep{tz-^), 

fe=i 

and say that a meromorphic hmction f(z) is p-elliptic if f{pz) ~ f{z). The simplest 
nonconstant p-elliptic function (of the second order) has the form 6p{az, bz)/6p{cz, dz), 
where ab = cd. 

The standard elliptic gamma function, depending on two complex bases p and q 
lying in the unit disc, \p\, \q\ < 1, has the form: 



rp,,(^) = n 

j,fc=0 



1 - z- V+ig'=+i 
1 — zp^q^ 



where z G C*. It obeys the properties Tp q{z) = Tg_p{z), 

^P.qiqz) = Sp{z)^P,q{z), ^p.qiPz) = Og{z)Tp,g{z) , 

and has zeros at z = p^>aq^>a and poles at z = p^<og^<o. The reflection formula 
has the form Tp^g(a)Tp,q(b) — 1, ab = pq; we set also 



Tp^q{ai, . . . , a„j) := ]^ Tp^q{ak), Tp^q{tz'^'^) := Tp^q{tz)Tp^q{tz 
fc=i 

rp,,(tzi'=^Z^^) = Tp,q{tZlZ2)Tp^q{tZlZ2^)Tp^q{tz^^ Z2)Tp^q{tZi^ Z^-^). 

2. The elliptic hypergeometric equation 

The following elliptic analogue of the Gauss hypergeometric function was intro- 
duced in [HITS] 

rp.glijZ^^) dz 



(1) V{t;p,q) = K . ^ ±2^ ^ 

where k — (p;p)oo{q',q)cx}/'2 and T is the positively oriented unit circle. The pa- 
rameters tj are restricted by the balancing condition Ilj^i — {PO)^ ^^d the 
inequalities \tj\ < I, j — 1, . . . ,8. The 1/- function can be meromorphically contin- 
ued to all tj G C* preserving the balancing condition. For tYts = pq (and other 
similar restrictions), it reduces to the elliptic beta integral [13] 

fo\ / 11^=1 rp,g(ijZ±i) dz Tl -r ^ \ 
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The addition formula for elliptic theta functions written in the form 

yields the following contiguity relation, via a corresponding relation for the inte- 
grands: 



hViqh) ^ t2V{qt2) ^ hV{qh, ^ ^ 
9p{tit2^ .tit"^^) Op{t2tf^,t2tf^) 9p{t3tf^ ,t3t2^) 

where V{qtj) denotes the (7)-function with the parameter tj replaced by 

qtj (with the balancing condition being Oj=i ^ P^l)- Bailey- type symmetry 
transformations for the ^-function [lOl E] give to ([3]) different forms, including in 
particular: 



UU(^pitltj/q)V{h/q) UU(^p{t2t,/q)V{t2/q) , 11 "=4 fe^,/?) ^^(^3/?) 



tldpih/tl^ts/tl) t20p{ti/t2,t3/t2) i3^p(il/i3, ^2/^3) 

where Oj^i — P^q'^- In combination with this yields the elliptic hypergeo- 
metric equation [15] : 

(5) A{ti,t2, . . . ,ts,q;p)(u{qti,q-'^t2;q,p) - U{t;q,p)^ 

+A{t2,ti, . . . ,ts,q;p)(u{q-Hi,qt2,;q,p) - U{t;q,p)'^ + U{t-q,p) = 0, 
where we have denoted 

8 



Op{tijt2,t2/qti,tit2lq) ^J-^ Opitstk) 



and 



U{t;q,p) := 



^p{ti/qh ,hti,h/t i) -p-p 9p{t2tk/q) 
2/qti,tit2 

V{t;q,p) 



Y[k=l Ip,i?(^fe^3 ) 

The potential A(ti, . . . ,ts,q;p) is a p-elliptic function of parameters fi, . . . , tg, one 
of which should be counted as a dependent variable through the balancing condition. 

We set ti = (pq)'^/t2 ■ ■ - tg, and perform the shift t2 pt2 (so that ti ti/p). 
Since the function A is p-elliptic in all parameters, we have A{p~^ti,pt2, . . .) = 
A{ti,t2, ■ ■ ■)■ The function U{p~^ti,pt2) defines therefore an independent solution 
of the elliptic hypergeometric equation. Let us compute the Casoratian of these 
two solutions (i.e., a discrete version of the Wronskian). For this, we multiply the 
above equation by U{p~^ti,pt2), the equation 

A{h,t2, . . .h,q]p){u{p'^qti,pq-h2) -U{p^hi,pt2)) 

+ A{t2MM, ■ ■ ■ ,q;p)(u{p'^q'hi,pqt2) ~ u{p-hi,pt2)) +u{p-hi,pt2)^Q 

by U{ti,t2), subtract them and obtain 

(7) A{h,t2, ...ts, q\p)D{p-Huq-H2) = A{t2, tuts,---, q;p)D{p-\-Hut2), 
where 

D{ti,t2) = U{qpti,t2)U{ti,pqt2) - U{qti,pt2)U{pti,qt2) 
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is the needed Casoratian. It is symmetric in p and q, which is an important property. 
The expression 

^^^^ ^ V{pqti,t2)V{ti,pqt2)-ti\^ V{qti,pt2)V{ph,qt2) 

I\l=i^P,g{tktf^^Pqtktf^) 

can obviously be interpreted as the determinant of a particular 2x2 matrix whose 
elements are expressed via the ^-function. 

Since t\ is a dependent variable, this is actually a first order difference equation 
in t2- After scaling ti pti, t2 — > 9^2 (so that ti = P?/ 11^=2 '^^ obtain the 
following equation for f{t2) ■= ^2): 

A{pti,qt2,h, . . ■,q;p) ,,, . 
A{qt2,pti,t3,...,q;p) 



h Op{ti/qH2,h/qt3,t^Hf') -pj- 9p{t2tk) . 
'qt2 9p{t2/tut2/ts,q-%Hi') O^ihtk/qy^ 



which yields 



D{ti,t2) - C(t2) , , , , T. /, ,±1N ' 



k=l 

where C{qt2) — C{t2). We can repeat the whole consideration with permuted p 
and q and obtain C{pt2) — C{t2)- This means (for incommensurate p and q) that 
C does not depend on t2, but it may depend on other parameters t3, . . . jtg. To 
compute C, we apply the residue calculus. For this we take the parameter from 
inside the unit circle to its outside and impose the constraints l^al > 1 > Igta], Ipts]. 
Then we deform the contour of integration T entering the definition of V{t) to the 
contour Tdef deformed in such a way that no poles are crossed during such a change 
of ts. The Cauchy theorem leads to 

v{t) := v,,^,it) = v,it) + n.tyr..fe^3^^) _ 

We take then the limit t2 l/ts and find the value of C through the limit for 
ratios of the left and right-hand sides of the above equality 

^ V{pqti,t2)V{ti,pqt2) - tr%"' V{qti,pt2)V{pti,qt2) ^ . ,. . 

For t2 1/^3, the function V{ti,pqt2) reduces to the elliptic beta integral, the 
residues of V{pqti,t2) blow up with VT{pqti,t2) remaining finite, the residues of 
V{pti, qt2) and V{qti,pt2) remain finite as well as the functions VT{pti, qt2) and 
VT{qti,pt2)- As a result, only the first term of our Casoratian survives and yields 

^_ n3<j<fc<8 ^P.9(^i^'=) 
^p,g{tl ^2 ) 

We obtain thus the formula 

(8) V(pqtut2)V{tuPqt2) - t^^t^^V{qh,pt2)V{pU,qt2) = ^'-^^^f J^'fi^"^ ■ 



^P,q\'-l H 
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The described solution of the elliptic hypergeometric equation is defined for 
l^l < 1, although equation ^ itself does not demand such a condition. It can be 
verified that 

A ( ■ • ■ ; ^A{ti,...,ts,q ^■,p) . 

The transformation tj p^^"^ /tj, j = 1, . . . , 8, maps therefore the elliptic hyperge- 
ometric equation to itself with the base change q —> q^^. Equivalently, the same 
inversion g — > 1/g occurs after the transformation tj p""^ /tj with integer such 
that X]j=i =4. As a result, we obtain the following solution of the elliptic 
hypergeometric equation in the regime |g| > 1 

rj,, , V{p^'Vt,,...,p^l^/t,-q-\p) 

\{k=l^P.q-^'yPl^kt3M/tk) 

As to the unit circle case \q\ — 1, the corresponding solution of equation (O can 
be obtained with the help of the modified elliptic gamma function or the modular 
transformation TS]. 

3. A CHARACTERIZATION THEOREM FOR THE V-FUNCTION 

We would like now to present contiguous relations for the ^-function (and, so, 
equation ([5])) in a 2 x 2 matrix form. For that we introduce the function 

y(ti, . . . ,t8;p, q) 



W{ti,...M:z) 



where z is some auxiliary variable. Replacing parameters ii,2,3 by ^4,7,8 and the 
function by in ([3]), we obtain after shifting t^ qts 

(10) W{qt3,qt7) = a{t; z)W{qt3, qt^) + (3{t- z)W{qt^, qts), 

where W{qtj, qtk) means the W{t; 2:)-function with respective parameters tj and tk 
replaced by qtj and qtk, and 

^ eAuz^\trtt^) ^ 0Auz^\t7ti^) 

' ep{t7Z±\tdt) dp{t7Z±\tstt') 

are p-elliptic functions of all variables (including z). 

Replacing now ^1,2,3 by (7^1,2,3 in O, and then permuting ti and t2 with t^ and 
<8, we obtain 

(11) W{qt3, qti) = -fit; z)W{qt3, qtg) + S{t; z)W{qt4, qtg), 
where 

, , _ 9p{tsZ^^,t3tg'^) -pr 6p{titj) 



'p(i4z±i,t3tr'),=i;^2 5,6.7 

(12) 5(i;^)^M^ff^M^ TT 



9p{t3z±\t4t^^) ^^^\\^^^.^ep{t^t,) 

are, again, p-elliptic functions of the parameters. Eliminating W{qt3, qt^) from (flQ] 
and lU), we obtain the relation 

(13) W{qt:i,qt7) = (a(i; z)7(i; z) + /3(t; z))W{qt3, qts) + a{t; z)Sit; z)W{qtA, qts). 
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We define now the matrices 

(14) M(tiM-MM-MM.tiM) 

(15) A{ti,t2'-,t^-,ti]tr,,tQ,tT,ts) 

where 



W{pti,qh) W{pt2,qh) 
W{ph,qU) W{pt2,qU) 

All Ai2 
A21 A22 



t5,...,ts; z;p,q) = a{pti,q^'^ts)j{pti,q~'^ts) + Piph, q~^ts), 
t5, ■ ■ ■,t8;z;P,q) = a{pti,q^'^ts)S{pti,q^^ts), 

t5, ■ ■ ■,t8;z;p,q) = ^12(^1, t2;i4,t3;i5, • ■ ■,t8;z;p,q), 
is, . . .,ts;z;p,q) = Aii(ti, t2; ^4, ^3; ^5, • ■ ■,h]z]p,q) 



^12(^1, ^2; 
^21 (il, ^2; ^3i 
^22(^1, ^2; ^3j ^4; 

are p-elliptic functions. In particular, 

Aij{p^^ti,pt2\h,U]h, ■ ■ .,ts;z;p,q) = A^j{ti,t2]h,U]h, . . . ,ts; z;p,q). 

After replacements t-j t-jx, tg, igx"^, equations (|13p and its partner obtained 

after permuting ^3 with are rewritten as a linear first order matrix g-difference 
equation 

(16) M{qx) = A{x)M{x), 

where we indicate only x-dependence. 

After the permutations p ^ q,ti is , ^2 ^4 , the matrix M is transformed to 
its transpose . Equation (fT6|) gets therefore transformed to 

(17) M{px) = M{x)B{x), 
where 

f B B 

(18) B{x) := B{ti,t2\h,ti\t^,tQ,^x,hx^^) -.^ [ " 

\ i>21 i>22 

Bii = ^11(^3,^4; ti,t2\h,h, irx, hx^'^]z] q,p), 
B12 = A2i{t3,t4;ti,t2;h,te,t7X,tsx~^;z;q,p), 
B21 = Ai2{U,U;ti,t2]h,U,t-jx,hx''^;z;q,p), 
B22 = A22{t3,t4;ti,t2;t5,t6,t7X,tsx~^; z;q,p). 



In principle, from the existence of first equation (jT6|) . it follows that there exists 
some p-difference equation of the form (|17|) with g-elliptic coefficients [6]. Indeed, 
we can simply take B{x) :— M{x)~^M{px) and see that 

B{qx) = M{xy^A{xy^A{px)M{px) = M{xy^M{px) ^ B{x). 

However, this i?-matrix is not unique. We let g{x) denote a matrix satisfying 
g{qx) = g{x). Equation (fT6|) does not change after the replacement M Mg, but 
the i?-niatrix gets changed to _B — s- g(x)~^B g{px) showing a functional freedom in 
the definition of this matrix. 

We suppose now, that M' is another meromorphic solution of equations pB]) and 
(flT)) . The matrix N = M'M~^ satisfies then the difference equations 

N{px) = N{x), N{qx) = A{x)N{x)A{x)~'^ . 

The first equation states that iV has p-elliptic entries, and the second one is a q- 
difference equation with the p-elliptic coefficients. The normalization chosen above 
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works for q- and p-shifts of all parameters, not just those of and ts, because of 
the permutational symmetry. 

Theorem 1. For incommensurate p and q, N is a constant multiple of the identity. 

The proof of this Theorem is given in a more general context in the last section. 
This statement simply means that the non-trivial p-elliptic functions (the matrix 
elements of N) cannot satisfy this g-difference equation with p-elliptic coefficients. 

With the help of three pairwise incommensurate quasiperiods Wi_2,3 and the 
parametrization 

q = e ^2 , p — e , 

we can convert q- and p-difference equations into linear finite difference equations. 
For this it is sufficient to pass to parameters gj introduced as tj = e^'^'^^/'^^, j = 
1, . . . ,8, and x = e^'^™/'^^^ function of u, the F-function is thus characterized 
as a unique (up to a constant independent of u) solution of 5 linear fininite difference 
equations: two difference equation in (|16p working with the shifts u u + uji, their 
partners in (jl7p working with the shifts u ~^ u+tos, and the condition of periodicity 
under the shifts u — > u -\- uj2, equivalent to the analiticity condition in x G C*. In 
order to characterize the proper function V{t; q,p) as a function of ii, . . . , tg up to 
a constant, we simply need to adjoin the permutational symmetry group Ss- 



4. Determinant representation of elliptic Dixon integrals 

The elliptic Dixon integrals (a.k.a. Type I integrals with BC'n symmetry) have 
the form: 

r 1 " T-r2n+2m+4 p ±1n 

4 (h,...,t2n+2rn+i}-^nl [[ .^il^il^l r (z^^) 2nJ^z ' 

^ l<i<J<n ^1 I j=l '^P.Qy^J I 



where \tj \ < 1, 



\m+l _ _ [PtPJooKIjDc 



By convention, when n — 0, /q™'' '■— 1; when n = 1, the resulting univariate integral 
is a higher-order version of the elliptic beta integral. 

The following transformation identity has been proved in [lOj . 

Theorem 2. The integrals ll^^ satisfy the relation 
(19) 

■ ■ ■ ,t2n+2m+A) ^ TT ^P,q{^rt s) [ — 

l<r<s<2n+2m+4 ^ 

In particular, when m = 0, the integral on the right-hand side is 0-dimensional, 
and one obtains an explicit evaluation of the left-hand side. 

One of the objectives of the present work is to give a more elementary proof of 
this transformation, by showing that both sides satisfy the same family of difference 
equations and initial conditions. We will also obtain a new proof of the special case 
m — 0. 
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The key idea is that can be written as a determinant of integrals of the 
form /^"+™ j TffQ -will thus be able to use difference equations for the univariate 
integrals to deduce difference equations for 7,^™'' . 

Due to the reflection identity Tp,q{z)Tp,q{pq/ z) = 1, we can write the /i™^ inte- 
grand's cross factors as 

n FTik^- n (^r^^p(-.-f ) -r^^.(-.-f )) ■ 

l<i<j<n '^P'ly'^i ) l<i<j<n 

The point, then is that the antisymmetric factor Y[i<i<j<n -^r^^pi-^i-^f^) '^^^ 
written as a determinant; more precisely, we have the following elliptic analogue of 
the Cauchy determinant: 

(which, as observed in can be obtained from the usual Cauchy determinant via 
a suitable substitution), using which we can bring our integral to the form 

n+2m+4 / 



r " /T-r2n+2m+4p /, ±1n „ 



j=i \ -Lp.gl^j ) k=i iZj^ 

X det (zj) det V* (^i ) , 



where (j)i{zj) = ai/9p{aiZj ), ipi{zj) ~ bi/9q{biZj ). But then, using the Heine 
identity: 



"T [ .J^^t Mzj)^det iptizj) TT = det / (f>i{z)i:j{z)dn{z), 

l<i<n 

an write 

-^n (^1 ) • ■ • I i2n+2m+4) = JJ^ 



we can write 



l<^,<n y A r,,,(z±2) li JJ. 2vry^. 

where k = (p;p)oo('Z; 9)oo/2. If we choose a.; = ti, bi — tn+i, 1 < i < fi, then the 
entries of the above determinant are of the form ^\ To be precise, the ij 

entry is 

Tq{ti) ^Tp{tn+j) ^ \qti, . . . ,qtn,ptn+l, ■ ■ ■ ,pt2n,t2n+l, ■ ■ ■ ,t2n+2m+4), 

where Tq{tk) represents the g-shift operator, Tq{tk)f{tk) — f{qtk)- We can also let 
the sequences a, b overlap, at the cost of some slight extra complication. For in- 
stance, for n — 2^m — and ai — bi = U, i — 1,2, we obtain on the right-hand side 
the Casoratian of F-functions which yields I^\ti, . . . , is) = ni<j<fc<8 ^p,q(tjtk)- 

This result hints that the integral /n^"* is computable in the closed form for arbitrary 
71 > 2. 
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Note that when m = —1, the balancing condition reads ti ■ ■ ■t2n+4 = 1, and 
thus at least one of the entries of the determinant cannot use the unit circle as its 
contour. This is not, however, a serious issue, since there always exists some valid 
choice of common contour, and the Heine identity works regardless. 

Now, the fact that the integral vanishes for m = — 1 (a fact visible from the 
explicit formula for m = 0) implies that the rows of the above matrix must be 
linearly dependent. Conversely, if we can find a linear dependence between integrals 

Tq{ti) ^Tp{tn+j) ^Ii^ '^\qti, . . . ,qtn,ptn+l, ■ ■ ■ ,pt2mt2n+l,t2n+2), 

which is independent of j (i.e., is p-elliptic in . . . ,<2n), that will imply van- 
ishing for m = —1, which as shown in [4] allows one to compute the integral for 
m = 0. This leads us to examine recurrence relations for the univariate integrals. 

In general, not only are the entries of the above matrix integrals of the form 
j{n+m 1)^ ^^^^ k X k minors themselves are proportional to integrals of 

the form /^"^"^ The recurrence for univariate integrals implicit in the vanishing 
of 4"^^ gives rise to recurrences for higher /i"*^ in the following way. 

Lemma 3. Let M be a n x k matrix, and suppose the vector v satisfies vM = 0. 
Then the k x k minors of M satisfy the (ri — k + l)-term relation 

y V, det (Mil,) = 0. 

le{l k~l.i},l'£{l....,k} 

k<i<n 

Proof. This certainly holds, by linearity, if we were to extend the sum down to i = 1, 
but the additional terms all have repeated rows in the minors, thus vanish. □ 

There are two main sources of recurrences for hypergeometric integrals. The 
first is recurrences of the integrands themselves (so long as the contour conditions 
can be satisfied by a common contour, that is). 

Theorem 4. The integral In^^^ satisfies the {n + 2) -term recurrence 



(20) 

'\ ^^, T,{t,)l!rHh, . . . , W2„.+4) = 0, n *. = (m)>. 



2n+2m+4 



/ y i-rn+2 

l<i<n+2 



Proof. If we divide out by common factors of the integrand, this reduces to the 
relation (see Lemma A.l in [3] or Corollary 2.3 in [3]): 

l<i<n-|-2 lli=l, '^Pv'-«''j ) l<j<n 

□ 

This approach alone is insufficient to get the full system of difference equations; 
for one thing, it is completely independent of the balancing condition. As in [9], 
the key is to multiply the integrand by functions related to the difference operators 
of [I0|. 

With this in mind, we consider the following function: 
g^")(2;;ii, . . . ,t„+2;t;i, . . ■,Vm+4) 

^ Ul<^<m+i^pM ni<.,<rn+4^p(^»/^) 

l<i<m+2^piiPl/^i)^) Z ^9p{z ^) ni<i<m-|-2 ^p((P9/^i)/^) 
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By inspection, this is invariant under the change z^Xj z. In addition, if ni<i<rn+2 
ni<i<m+4^» ~ (P?)'""''^?) then both terms are (meromorphic) p-theta functions 
with the same multipher, and thus 

g'^"'\pz;tl, . . .,tm+2;Vl, . . . ,Uto+4) = pz'^g'^"'\z;tl, . . .,tm+2;Vl, . . . , 

Moreover, the apparent poles at z = ±1, ±^/p must by symmetry have even order, 
and thus g^"^^ must in fact be holomorphic at those points. Thus g^"^^ has only 
simple poles at the points {pq/ti)^^p^; it follows that 



g^"'\z;ti, . . .,tm+2;Vi, . . .,Vm+A) = ^ 



l<i<m+2 P^^^ ' ^' 

for suitable coefficients Ui which can be computed in the usual way: clear the 
denominator and set z = tilq to find 

The relevance of this function for our purposes is that it integrates to against 
the 7^"*^ density. More precisely, we have 



/ Am)(^.+ f. .f. f. ■^ Ill<r<2m+6^P,qitrZ^^) dz 

/ 9 {Z,ti, . . . ,tm+2,tm+3, ■ ■ ■ ,t2m+6) f; , +n. r 7=^—^^ 

i|z|=l >^p,q{z ) 2-K\/-\z 

SO long as Wx<r<1m+%ii = (w)""^^? and \tl\/q, . ..,\tm+2\/q, \tm+3\, |i2m+6| < 

1. Indeed, by symmetry, both terms of g*^™) have the same integral; if we restrict 
to the first term (thus gaining a factor of 2), perform the change of variables z i-^ 
q~^^'^/z, and move the contour back to the unit circle (which crosses over no poles), 
we obtain 

But now the integrand is antisymmetric with respect to ^ i— > and therefore the 
integral vanishes. 

Using the partial fraction decomposition of we thus obtain a new recur- 

rence. 

Lemma 5. If h ■ ■ ■ t2m+6 = {pq)"^~^^q, then we have the following {m + 2)-term 
recurrence for ' : 

nm+3<»<2TO+6^pft'^fc/g) rr. u s-lj{m.),. . n „ 

(21) 2^ — Ff n I. -. Tqitk) !{ '{h,...,t2m+6) = 0- 

l<k<m+2 '^fc lll<i<m+2;i7^fc ^pKH/J-k) 

Applying the operator Tp{tm+2+i)~^ to this equality, we obtain 

/r,o^ n Ilm+3<i<2m+6^p(^i'^k/q) 

(22) 2^ vkMki = 0, vk = — rrriTY' 

l<fe<m+2 lll<i<m+2;ij!:k''pV^i/'^k) 
Mkl = Tq{tk)~^Tp{tm+2+l)~^ Ii^\tl, . . . ,i2TO+6)- 

Corollary 6. Ift\-- ■t2n+2 = 1, then In = 0. 
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Indeed, as shown above the integral In is proportional to the determinant 



of the matrix AI in (|22p with m = n — 2 and scaled parameters. However, the 
vector V = (ui, . . . ,ti„) belongs to the kernel of M, vM = 0, and, so, detM = 0. 
For n = 2, such a result follows also from ([8]) and the elliptic beta integral. This 
statement proves the vanishing hypothesis of [1| , which was needed there for a proof 
of the evaluation formula for the elliptic Dixon integral. 

Applying Lemma [3] to relation (|22p . we obtain equality X^fc^f '^fe^fc — 0; where 
dk = det{Miii) with I e {1, . . . , n — 1, A;}, I' e {!,..., n}. The minors dk are pro- 
portional to the integrals In'^^'"'^^\ti/q, . . . , i„_i/g, . . . , ifc/g, . . . , tm+2,tm+3/p, 
. . . , tm+3+n/p, tm+i+n, ■ ■ ■ , t2m+&) ■ Substituting corresponding explicit expressions, 
multiplying parameters ii,...,t„_i by q and t^+a, ■ • ■ , *m+3+n by p, changing 
TO ^ m + n — 1 , and permuting parameters appropriately, we obtain a recurrence 
for general -integrals. 

Theorem 7. If ti ■ ■ ■ t2m+2n+4 — t/ien we have the following (m + 2)- 

term recurrence for /i™'' ; 

/'OQ^ nm+3<t<2n+2m+4^p(^»^fc/g) Tn n „ 

(^■^j 2^ ~Tn n 1+ 1+ \ J„ (El, ■ • ■ ,l2«+2m+4j = U. 

l<fe<m+2 lll<i<m+2;i5^fe "pV-il'-k) 

Corollary 8. [MIIS] //ti • • •t2n+4 = pg, </ien 

l<i<j<2n+i 

Proof. Indeed, both sides satisfy the same 2-term recurrence, and thus their ratio 
is invariant under Tq{ti)~^Tq{tj), and similarly for p-shifts. It follows that their 
ratio is independent of ti, . . . , i2n+4- To determine the remaining factor, we may 
consider the limit of the ratio as ^271+3^271+4 ^ 1, and proceed by induction. □ 

Corollary 9. Ifti-- ■t2n+2m+A = PQ, then we have the {"^^^ -dimensional deter- 
minant 

det ( TT Tp{tr) TT Tg{ts+n+m)li"'\ti, . . . , t2„+2m+4) ) 

i?,SC{l,2,...,n+m};lfl| = |S|=m I \\ " / 

\reR sGS / 



[ JJ^ tj(^p{titj )tn+m+jOp{tn+m+itn+m+j)\ 

l<2<j<n+m 

]^ ^p,q{tit 



X 



l<i<j<2n+2m+4 



Proof. Using the determinantal representation of I^™-* , we can express the (R, S) 
entry of the above determinant as 

tjdp{titj ) tn+m+j(^p{tn+m+itn+m+j . 

i<j i<j 

^ ."^tX n n 7^«(Wm+.)4"^""'Hil,---,W2m+4)), 

j^S^ l<r<n+m l<s<n+m 
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where i?^ = {1,2, . . . ,n + ni} \ R. The first two factors can be pulled out of the 
determinant, as they are independent of the column or row as appropriate; this 
gives an overall factor 

n n (tMutf))-' = n n itjOpiutf))-' 

RC{l,2,...,7i+m} iJeR" l<i<j<n+m i?C{l,2,...,n+m} 

\R.\=m i<j \R\=m- i,j(^R 

= n wt^tf^y^"^"-'^'^- 

l<z<j<n+m 

We are thus left with computing a determinant of determinants. These are all 
minors of a fixed matrix, and thus our ("^™) x matrix is the rt-th exterior 

power of the n + m x n + m matrix with ij entry 

n+m+r )Il {tl, ■ ■ ■ , t2n+2Tn+4) , 

l<r<n+m 

i.e., the matrix giving the action of the original matrix on the rt-th exterior power 
of its natural module. In general, the n-th exterior power of the n + m x n + m 
matrix M has determinant 

det(Af)(""-i"'), 

which can be seen easily by reduction to the diagonal case (the naturality of the 
construction implies invariance under conjugation). In our case, the n + mx n + m 
matrix has determinant 



{tl,...,t2n+2rn+4), 

l<i<j <n-\-rn 



which implies the desired result. □ 

Remark 1. When m = 1, this is essentially an elliptic version of Varchenko's 
determinant of univariate hypergeometric integrals [El [12]. Indeed, Varchenko's 
determinant is equivalent to Dixon's integral evaluation, which is a limit of the In^ 
evaluation formula 11, Thm. 7.2]. Convergence of the matrix itself is somewhat 
more subtle, as the domain of integration tends to pass through algebraic singular- 
ities of the integrand in the limit, thus giving rise to somewhat tricky phase issues. 
One does find, however, that 

lim '— - - 

q'^i ai,q°'^ /ai, g"^ aa , q"^ /as ; p, q) 

|^?p(aia^^)r-"-"^^J^^i^(2-(p;p)L) 
r(ai)r(a2) 



(24) 



Jze[ai.a2] l<r<n ZTT^/-iZ 



where ai, a2, bi,. . . ,bn are on the unit circle with positive imaginary part and 
3fi(ai) > 3fi(a2), the exponents satisfy the convergence conditions ^{af) > > 
5)?(/3~), and one has the balancing condition ai+ a2 — J2i Pii with ai := al + , 
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Pi := /\ + f3j^ . The domain of integration is the counterclockwise arc from ai to 
a2- The proof is as in Theorem 7.2 of [11]; note also that by Lemma 7.1, op. cit., 
one has 

OpibrZ^y- = exp{-(3r{log{-brZ)+\og{-br/z))/2)\0p{brZ^^)f^. 

The resulting factor is in particular locally constant in z on the upper semicircle, 
changing only when z passes over one of the singularities bi. Also, as noted in ll], 
the change of variables x — —9p{z)'^/9p{—z)^ turns this into a higher-order ordinary 
beta integral, precisely as appears in Varchenko's determinant identity. Our result 
can also be viewed as a generalization of the main result of the work of Aomoto and 
Ito [3] , which corresponds to a further degeneration of the trigonometric integral 
of [m Thm. 5.6]; the Jackson integral is obtained as the sum of residues obtained 
upon shrinking the contour to 0. 

The integrals /n™'' thus provide solutions of the system of recurrence relations 
([20| and ([23|l and their partners obtained after permutations of parameters. How- 
ever, these recurrence relations do not imply the constraint l^l < 1 which is needed 
for the definition of The construction of solutions of the elliptic hypergeo- 

metric equation with l^l > 1 ([9]) extends to arbitrary values of n and m. Modulo 
some ellipticity factor, the g-shift equations satisfied by 

l<i<2n+2m+4 

are invariant under the transformation ti i-^ p""^ /ti^ q i— > 1/g, where oi, . . . , a2n+2m+4 
is any sequence of integers or half-integers (but not mixed) such that ^2n+2m+4 ^ _ _ 
2m -f 2. That is, for both recurrence relations, the rescaled coefficient of each term 
gets multiplied by the same quantity under such a transformation. This observa- 
tion provides us with the solutions of those recurrences for \q\ > 1. Solutions with 
|g| = 1 are obtained (as in the n = m = 1 case) by using the modified elliptic 
gamma-function which we do not consider for brevity. 

5. A PROOF OF THE TRANSFORMATION FORMULA 

The two recurrences we have given for In^^ are also recurrences for the right-hand 
side of relation (I19p : the transformation simply swaps the two kinds of recurrence. 
The proof of Theorem [5] using these difference equations boils down to showing that 
/i™'' is the unique meromorphic solution of the full system of p- and g-difference 
equations, up to an overall constant (which is then easy to obtain via a limit). As a 
subcase, this proves Theorem[l]as well. It is not particularly elementary; it involves 
some difference/differential Galois theory [5]. We start by giving the key auxiliary 
statement needed for us. 

One technical issue that arises in our application of difference Galois theory is 
that the literature, and many of the main statements, require that the constant 
field be algebraically closed. Since we are dealing in our case with a family of 
difference equations, this constraint is too strict. However, we can evade this issue 
via a suitable base change, as follows. Let (fc,T) be a difference field (i.e., r is an 
automorphism of the field k) of characteristic such that the field k'^ of constants 
is algebraically closed in k. (In other words, the only finite orbits of r on fc are of 
length 1.) There is then a canonical extension of r to the field I := k (E)k^ k'^ , such 
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that the algebraic closure k'^ is the new constant field. We thus obtain a well-defined 
notion of Galois group over I. 

Lemma 10. Let A e GLn{k), and let W be the space of vectors w G k'^ such that 

t(w) = Aw. Let G be the Galois group of this difference equation (viewed as an 
equation with coefficients in I), with associated representation V. Then 

dimfex(VF) = dim((y<^). 

Proof. The absolute (ordinary) Galois group Gal(k^) acts naturally on I, and 
commutes with the action of r; in particular, we have a natural isomorphism 
Gal{k'^) = Gal{l/k), and the latter is well-defined. In particular, if wc extend 
coeflScients of the difference equation to /, the resulting vector space is stable under 
Gal{l/k), and thus admits a basis of fc-rational vectors. In other words, the dimen- 
sion dim/jT (W) is unchanged under this coefficient extension, and wc may therefore 
assume k = I, or in other words that the field of constants in k is algebraically 
closed. 

Now, let the fc-algcbra i? be a Picard-Vessiot ring of the difference equation. 
Then, by definition, the set of solutions w £ R" of the difference equation is an 
n-dimensional vector space over k'^ with a G action equivalent to the representation 
V. A solution has coefficients in k iff it is invariant under the action of G, and thus 
the result follows. □ 

Remark 2. A similar lemma holds for the differential case (replacing GL„ by its 
Lie algebra); note that in that case, the constant field is automatically algebraically 
closed in the coefficient field. 

Remark 3. More generally, if p is any rational representation of GLn{k), one can 
apply the lemma to the difference equation t(?«) = p{A)w, using the fact that the 
new Galois group is simply p{G). 

Applying remark 3 to the adjoint representation gives the following: 

Corollary 11. Let A.G be as above, and let W be the space of matrices M e 

End{k") such that t{M) = AMAr^. Then 

dimfe.(VK) = dimiiEndaiV)). 

In particular, if G is irreducible, then W is 1- dimensional, and (since I & W ) W 
consists of scalar matrices. 

That EndoiV) = I when V is irreducible is a standard fact (Schur's lemma) of 

representation theory. Any cndomorphism which is not a multiple of the identity has 
at least one proper eigenspace (since I is algebraically closed), and each eigenspace 
is invariant, making the representation reducible. 

We can pass now to the proof of our transformation formula itself. In this case 
the field k coincides with the field of p-elliptic functions in one of the parameters 
X (actually, in all parameters), and r is the (/-shift operator, T{f{x)) = f{qx). 
Therefore k'^ is a field of constants independent on x (the only simultaneously p- 
and g-elliptic set of x-functions). In order to apply the above theory, we need 
to show three things: first, that our two recurrences can be combined to give a 
difference equation in matrix form as above, second, that the coefficients of this 
difference equation can be made elliptic upon suitable renormalization, and third, 
that the resulting elliptic difference equation has irreducible Galois group. 
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Proof. To sec that the recurrences combine to give a matrix difference equation of 
order ("^'") , we need to show that one can use the recurrences to express 

(qti, qtm,tm+i, ■■■^qt t2m+2n+4:/ q', P, q) 

as a hnear combination of 

[]Tg(fi)4™)(tl, . . .,t2m+2n+4;P, q), 

ies 

where S ranges over m-element subsets of {1, 2, . . . , n + m}. Using the n + 2-term 
recurrence, we obtain a hnear dependence between the original integral, the integral 

I^Kltlj • • • ) qtm, tm+l, ■ ■■, t2m+2n+3, t2m+2n+i',P, 

and the n integrals 

m 1 tm+l, ■ ■ ■ ,^2 m+2n+3 7 ^2m+2n+4 /q\p, q) 

for m+1 < i < m + n. By symmetry, it suffices to consider the term with i = m+1. 
But then the m + 2 term recurrence gives a linear dependence between the m + 2 
integrals 

Tq{ti)~^ I^\qt\, . . . ,qtm+l, ■ ■ ■ ,t2m+2n+3,t2m+2n+4',P, Q) 

for ze{l,...,TO + l, 2m + 2n + 4} . Wc thus obtain a matrix difference equation of 
the form required. (One also notes that the corresponding matrix A is quite sparse; 
the entry corresponding to a pair S, T of m-subsets of {1, 2, . . . , n + m} is unless 
SnT >m-l.) 

Next, for ellipticity, we consider the renormalization 

' ' n [tl, ■ ■ ■ ,t2m,+2n+4:,P,q)- 



\Yll<r<2m+2n+i^P,qi^r'V^^) J 

This differs from the corresponding minor of the matrix for 7^"^+" by multipli- 
cation by a pair of factors, one of which is p-elliptic in all variables other than v, 
and the other of which is similarly g-elliptic. We thus find that the matrix A is 
essentially just the n-th exterior power of the corresponding matrix for 
up to a pair of diagonal matrices that combined have no effect on ellipticity. It will 
thus sufHce to show that the difference equation is elliptic when n = 1. 
Consider the m + 1 x m + 1 matrix 

X ll^'^ipqti, . . . ,pqtm+i,tm+2, ■ ■ .,t2m+A,x, 1/Tx;p,q), 

where T = ni<r<2m+4*^- "^^^^ symmetrical between p and q (being replaced 
by its transpose when p and q are swapped), so it suffices to consider its behavior 
under a q shift. One finds, in fact (using the two recurrences as described above), 
that 

M{qx) = A{x)M{x), 
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where 



6p{xv±\Txtf^) ^ 



9p{U{Tx)±\xv±^)) 



9p{Tx^)Ui 



<l<m+ 



,ep{tiTx) epitjx)ii 



m+2</<2m+4 



9p{l/titj) 



<l<m+l;l^j 

Op{ti/t,) 

These coefBcients are readily verified to be elliptic as required. 

It remains only to prove that the Galois group is (generically) irreducible for all 
TO, n. Now, the Galois group for n > 1 is the n-th exterior power of the group for 
n = 1, so it will suffice to show that the generic group for n = 1 contains SLm+i 
(all nonzero exterior powers of which are irreducible). We can proceed as in the 
proof of Theorem 3.3.3.1 in [1 . The point is that, by Andre's theory, the Galois 
group can only become smaller under specialization, including degeneration to a 
differential equation. We have already discussed the fact that the integral l}'"'' can 
be degenerated to a higher-order classical beta integral, and indeed one can obtain 
a basis of the corresponding differential equation in that way. It thus follows that 
the elliptic hypergeometric difference equation degenerates under that limit to the 
Jordan-Pochhammer differential equation. But this is known [TH] to have generic 
Galois group containing SL„i+i- □ 



Remark 4- Note, in particular, that the given formula for A indeed converges 
to the identity matrix in the Jordan-Pochhammer limit, as in particular T — > 
1, making the diagonal contribution converge to 1; the ofF-diagonal contribution 
vanishes since tjtj+m+i 1. One can moreover directly compute the limiting 
differential equation, and use the rigidity of the Jordan-Pochhammer equation to 
verify that the two equations are equivalent. 
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